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^ ' Random matrix models provide an interesting framework for modeling a num- 

ber of physical phenomena, with applications ranging from atomic physics to quantum 
gravity 2 . In recent years, non-hermitian random matrix models have become in- 
creasingly important in a number of quantum problems 3 ' 4 . A variety of methods have 
been devised to calculate with random matrix models. Most prominent perhaps are 
the Schwinger-Dyson approach 2 and the supersymmetric method 5 . In the case of 
Non-hermitian Random Matrix Models (NHRMM) some of the standard techniques 
fail or are awkward. 

In this talk we go over several new developments regarding the techniques 6 ' 7 for a 
large class of non-hermitian matrix models with unitary randomness (complex random 
numbers). In particular, we discuss 

(a) - A diagrammatic approach based on a 1/iV expansion 

(b) - A generalization of the addition theorem (R-transformation) 

(c) - A conformal transformation on the position of pertinent singularities 

(d) - A 'phase' analysis using appropriate partition functions 

(e) - A number of two-point functions and the issue of universality. 

Throughout, we will rely on two standard examples: a non-hermitian gaussian 
random matrix model (Ginibre-Girko ensemble 8 ), and a chiral gaussian random matrix 
model in the presence of a constant non-hermitian part 9 . The first ensemble being 
a text-book example will allow for a comparison of our methods to more conventional 
ones, the second ensemble will show the versatility of our approach to new problems with 
some emphasis on the physics issues. Further applications will be briefly mentioned. 



* Talk presented by MAN at the NATO Workshop "New Developments in Quantum Field Theory" , 
June 14-20, 1997, Zakopane, Poland. 



1 



— i 

z 

b 



in 



Figure 1. Large JV "Feynman" rules for "quark" and "gluon" propagators. 



DIAGRAMMATIC EXPANSION AND SPONTANEOUS BREAKDOWN 
OF HOLOMORPHY 

The fundamental problem in random matrix theories is to find the distribution 
of eigenvalues in the large N (size of the matrix M) limit. According to standard 
arguments, the eigenvalue distribution is easy to reconstruct from the discontinuities 
of the Green's function 

where averaging is done over the ensemble of N x N random matrices generated with 
probability 

P( M ) = -e~ NT ^ M \ (2) 
Z 

To illustrate our diagrammatic arguments let us first consider the well known case of a 
random hermitian ensemble with Gaussian distribution. 



Hermit ian diagrammatics 

We use the diagrammatic notation introduced by 10 , borrowing on the standard 
large N diagrammatics for QCD 11 . Consider the partition function 

Z =< det{z -H)>= J d*fjd^dHe- c e-% TlH2 (3) 

with a "quark" Lagrangian C 

£ = $ a (zl b a - H b a )i> b , (4) 

where H is a hermitian random matrix with Gaussian weight ( the width of the Gaussian 
we set to 1). We will refer to ^ as a "quark" and to if as a "gluon". The "Feynman 
graphs" following from (f|) allow only for the flow of "color" (no momentum), since (§) 
defines a field-theory in + dimensions. The names "quarks", "gluons", "color" etc. 
are used here in a figurative sense, without any connection to QCD. The "quark" and 
"gluon" propagators (double line notation) are shown in Fig. [I]. 

Introducing the irreducible self energy S, the Green's function reads 

In the large N limit the equation for the self energy S follows from resumming 
the rainbow diagrams of Fig. |2|. All other diagrams (non-planar and "quark" loops) 
are subleading in the large N limit. The consistency equation ( "Schwinger- Dyson" 
equation of Fig. §) reads 

E = G. (6) 
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Figure 2. Diagrammatic expansion of Green's function ([!]) for Gaussian randomness. 




Figure 3. Schwinger-Dyson equation. 

Equations (|5]) and @ give immediately G(z — G) = 1, so the normalizable solution for 
the Green's function reads 



G(z) = -(z - v^4) (7) 

which, via the discontinuity (cut) leads to Wigner's semicircle for the distribution of 
the eigenvalues for hermitian random matrices 

i/(A) = ^v 7 !^. (8) 



Non-hermitian diagrammatics 

If we were to use non-hermitian matrices in the resolvent (HI), then configuration 
by configuration, the resolvent displays poles that are scattered around z = in the 
complex z-plane. In the large N limit, the poles accumulate in general on finite surfaces 
(for unitary matrices on circles), over which the resolvent is no longer holomorphic. The 
(spontaneous) breaking of holomorphic symmetry follows from the large N limit. As a 
result dG/dz ^ on the nonholomorphic surface, with a finite eigenvalue distribution. 
In this section we will set up the diagrammatic rules for investigating non-hermitian 
random matrix models. In addition to the "quarks" we introduce "conjugate quarks", 
defined by the + dimensional Lagrangian 

C = ^(z- M)^ + 4>{z- M^)<p. (9) 

For hermitian matrices, "quarks" if) and "conjugate-quarks" decouple in the "ther- 
modynamical" limit (N — > oo). Their respective resolvents follow from (|9|) and do 
not 'talk' to each other. They are holomorphic (anti-holomorphic) functions modulo 
cuts on the real axis. For non-hermitian matrices, this is not the case in the large N 
limit. The spontaneous breaking of holomorphic symmetry in the large N limit may 
be probed in the z-plane by adding to (|9|) 

C B = X{^ + ^) (10) 

in the limit A — > 0. The combination £ + £b will be used below as the non-hermitian 
analog of the Lagrangian (|J). 

From (^|,|10]) we define the matrix- valued resolvent through 

(11) 
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Figure 4. All "gluonic" amplitudes for complex gaussian matrices. 

where the limit iV — ► oo is understood before A — > 0. The "quark" spectral density 
follows from Gauss law 13 , 

u(z, z) = -d- z G{z, z) = -\-d- z Tr N Q qq (12) 

which is the distribution of eigenvalues of M.. For hermitian Ai, (fl2|) is valued on the 
real axis. As A —>■ 0, the block-structure decouples, and we are left with the original 
resolvent. For z — > -MO, the latter is just a measurement of the real eigenvalue distri- 
bution, as shown before in the case of Gaussian hermitian ensemble. For non-hermitian 
J\A, as A — >• 0, the block structure does not decouple, leading to a nonholomorphic re- 
solvent for certain two-dimensional domains on the z-plane. Holomorphic separability 
of @ is spontaneously broken in the large limit. For more technical details we refer 
to the original work 7 . Similar construction has been proposed recently in 12 . 

Examples 

Consider first the Ginibre-Girko ensemble, i.e. the case of complex matrices with 
the measure 



< . . . >= J [dM] exp 
The "gluon" propagators read 



A tr(A^A^ f - rReMM) 



1-r 



(13) 



(\M ab \ 2 ) = 1 (M ab M ba ) = ^ (14) 

corresponding to hermitian (r = 1), anti-hermitian (r = — 1) or general complex (r = 0) 
matrix theory. 

From (|9|, [l0|) we note that there are two kinds of "quark" propagators (1/z for 
"quarks" ip and 1/z for "conjugate-quarks" <j), where both can be incoming and outgo- 
ing). The relevant "gluonic" amplitudes correspond now to Fig. ||a-§]d, where the (c,d) 
contribution corresponds to twisting the lines with a "penalty factor" r. 

The equation for the one particle irreducible (1PI) self-energy follows from Figs. f§- 
|5| in the form 

/El S 2 \ = J_ tr / G qq Gqq 
\ S 3 S 4 J N \ Gqq Gqq 

= til) °{i ;)■ ^ 

Here the trace is meant component- wise (block per block), and the argument of the 
trace is the dressed propagator. The operation o is not a matrix multiplication, but a 
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Figure 5. Self-energy equation for non-hermitian matrices. 

simple multiplication between the entries in the corresponding positions. Here tr/v is 
short for the trace on the N x N block-matrices. 

Each of the entries has a diagrammatical interpretation, in analogy to the hermi- 



tian case. For example, the equality for the upper left corners of matrices in fllq) is 
represented diagrammatically on Fig. [5[ The first graph on the r.h.s. in this figure 
does not influence the "quark-quark" interaction - it corresponds to the double line 
with a twist, therefore, as a non-planar one, is subleading. However, this twist could be 
compensated by the twisted part of the propagator coming from the second correlator 
(|HD, thereby explaining the factor r in the upper left corner of (|I5fl. Other entries in 
fll5|) follow from Fig. f| by inspection. 

The "quark" one-point function is now 

G(z,z) = — tr N Q m = (z - E 4 )/de£. (16) 
It follows that E 2 = E 3 , with 

det ■ Ex = t{z - E 4 ) (17) 

det ■ E 4 = t(z - Ei) (18) 

det-E 2 = E 2 -A, (19) 

where det — (z — ^i)(z — E 4 ) — (A — S 2 ) 2 . Substituting r = E 2 — A in the last relation 
in (^) yields the equation 

((*-£i)(z-£ 4 ) -r 2 )(r + A) =r. (20) 

For A = 0, the solution with r = is holomorphic while that with r ^ is nonholo- 
morphic. In the holomorphic case, Ei(z — Ei) = r, and the resolvent is simply 



CM = i^gBl (21) 

where the upper sign corresponds to the solution with the pertinent asymptotics. In 
the nonholomorphic case, G(z, z) = z — E 4 , with 

7 — T7 

G(z,z) = Y —^ (22) 

in agreement with 13 . The boundary between the holomorphic and nonholomorphic 
solution follows from the condition E 2 = imposed for the nonholomorphic solution, 
here this is equivalent to \G(z, z)\ 2 = \G(z)\ 2 = 1, that is 

X2 ■ y2 1 (23) 



l + r) 2 (1 



T 



,2 



which is an ellipse in the complex plane. Inside (p3|) the solution is nonholomorphic 
and outside it is holomorphic. The case investigated by Ginibre 8 follows for r = 0. It 
is pedagogical to compare our method of solving this problem to the one coming from 
interpreting the Ginibre ensemble as a two-matrix (7i ~ H\ + iH 2 ) model. 
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As a second example, we consider the Chiral Random Matrix model, which got 
recently some attention as a schematic model for spontaneous breakdown of the chiral 
symmetry. Here we consider for simplicity the Gaussian version of such model in the 
presence of a non-hermitian part, "chemical potential" /x, as suggested by Stephanov 9 . 
The non-hermitian character comes from the property of Dirac matrices in Euclidean 
space. The form of the determinant stems from the constant mode sector of the massless 
and chiral Dirac operator at finite chemical potential 14 . The corresponding partition 
function reads 

Z =< det(z- w-M) > (24) 

where 

7 = ilo = ( ? ~n ) and M = ( n I ■ ( 25 ) 



The only novel features come from the "chiral character" of the matrix Ai, i.e. the 
fact that it anticommutes with the 75 = diag(lAr, —In)- Due to this fact, the "gluon" 
propagator T> inherits the block structure which in the tensor notation (see Fig. 1) 
reads 

V = ^( 7+ ®7-+7-®7+) (26) 

with 7± = (ljy±7s)/2 and the bare "quark" propagator 1/z gets modified to (z— /ry) -1 . 
As a result, the 1PI self-energy equations in the planar approximation are given by 

(27) 





where T> is the "gluon" propagator (p6|), and Sj are diagonal 2N x 2iV matrices. In- 
verting in ( p7| ) with respect to the "quark-conjugate quark" indices gives, after some 
elementary algebra, two kinds of solutions: 

(i) A nonholomorphic solution (£2 = S3 7^ 0) ( "quark-quark" resolvent) 



G(z, z) = ^tr N g„ = Z-iy- \~P^ (2? 



where z = x + iy, a result first derived in 9 using different arguments. 

(ii) For S 2 = S3 = we recover the holomorphic solution 9 ' 15 , Si (z) = G(z)l, 
£4 = S|, with G(z) fulfilling the cubic Pastur-like equation 

G 3 {z) - 2zG 2 {z) + (z 2 + /i 2 + l)G(z) -z = 0. (29) 

Note that in the case of this example the standard techniques of multi-matrix 
models do not apply. 



ADDITION LAWS 

The concept of addition law for hermitian ensembles was introduced in the seminal 
work by Voiculescu 16 . In brief, Voiculescu proposed the additive transformation (R 
transformation), which linearizes the convolution of non-commutative matrices, alike 
the logarithm of the Fourier transformation for the convolution of arbitrary functions. 
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This method is an important shortcut to obtain the equations for the Green's functions 
for a sum of matrices, starting from the knowledge of the Green's functions of individual 
ensembles of matrices. This formalism was reinterpreted diagrammatically by Zee 17 , 
who introduced the concept of Blue's function. Let us consider the problem of finding 
the Green's function of a sum of two independent (free 16 ) random matrices M.\ and 
M.2-, provided we know the Green's functions of each of them separately. First, we note 
that the 1PI self-energy E can be always expressed as a function of G itself and not 
of z as usually done in the textbooks. For the Gaussian randomness, E#(G) = G (see 
(H)). Second, we note that the graphs contributing to the self-energy T*i +2 (G) split into 
two classes, belonging to £i(G) and ^(G), due to the independence of probabilities 
P(M\) and P(M 2 ) and large N (planar) limit. Therefore 

S 1+2 (G) = S 1 (G) + S 1 (G). (30) 



Note that such a formula is not true if the energies are expressed as functions of z. 
Voiculescu R transformation is nothing but R{u) = Ti[G(u)]. The addition (|30|) reads, 
for an arbitrary complex u, Ri + 2(u) = R\{u) + R 2 (u). The R operation forms and 
abelian group. The Blue's function, introduced by Zee 17 , is simply 

B(G) = S(G) +G' 1 . (31) 



Therefore, using the identity G(z) — (z — £) \ we see that the Blue's function is the 
functional inverse of the Green's function 



B[G(z)} = z (32) 
and the addition law for Blue's functions reads 

B 1+2 (z) = B 1 (z)+B 2 (z)--. (33) 

z 



The algorithm of addition is now surprisingly simple: Knowing G\ and G 2 , we find (p2|) 
B\ and B 2 . Then we find the sum Bi +2 using (p3|), and finally, get the answer G\ +2 , 
by reapplying (|32]) . Note that the method treats on equal footing the Gaussian and 
non-Gaussian ensembles, provided that the measures Pi and P 2 are independent (free). 

The naive extension of this algorithm fails completely for the non-hermitian ma- 
trices. It is not a priori puzzling - the underlying mathematical reason for ( p3|) is the 
holomorphy of the hermitian Green's function, not fulfilled for the case of NHRMM, 
as demonstrated in the previous section. However, since we managed to extend the 
diagrammatical analysis to the NHRMM, it is still possible to generalize the addition 
formula using the parallel diagrammatic reasoning like in the hermitian case. The 
generalization amounts to consider the matrix- valued Green's function (pTT]) . The gen- 
eralized Blue's function 6 ' 7 is now a matrix valued function of a 2 x 2 matrix variable 
defined by 



B(Q) = Z = f " . (34) 




where A will be eventually set to zero. This is equivalent to the definition in terms of 
the self-energy matrix 

B{g) = Y, + g- 1 (35) 
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where S is a 2 x 2 self energy matrix expressed as a function of a matrix valued Green's 
function. The same diagrammatic reasoning as before leads to the addition formula for 
the self-energies and consequently for the addition law for generalized Blue's functions 

z = B 1 (g) + B 2 (g)-g- 1 . (36) 

The power of the addition law for NHRMM stems from the fact that it treats Gaussian 
and non-Gaussian randomness on the same footing 18 . 

Example 

Let us consider for simplicity complex random Gaussian matrices (r = 0), which 
we rewrite as the sum Hi+iH 2 , with Hi, H 2 hermitian. The generalized Blue's function 
for hermitian Hi follows explicitly from 

^(JJH-OiH (37) 

a matrix analog of @. The generalized Blue's function for anti-hermitian iH 2 follows 
from 



where the entries reflect the antihermicity (set r = — 1 in (|I5)). It is a straightforward 



exercise to check that the matrix equation (R3) with the generalized Blue's functions 



£>i and B 2 , corresponding to (|37|) and (0), reproduces two types of solutions ([H]) and 
(p2|) as well as the equation for the boundary ( here the circle) separating them on the 
z plane. 



CONFORMAL MAPPINGS 

The existence of the nonholomorphic and as well holomorphic domains in the case 
of two solutions of NHRMM provides a powerful way to evaluate the supports for 
the level densities of NHRMM. The envelopes of these supports (supports form two- 
dimensional islands) can be derived very generally using a conformal transformation 
that maps the cuts of the hermitian ensemble onto the boundaries of its non-hermitian 
analog. 

Let us consider the case where a Gaussian random and hermitian matrix H is 
added to an arbitrary matrix M. The addition law says 

R h +m{u) = Rh{u) + R M {u) = u + R M {u). (39) 

where we have used explicitly that for Gaussian Rh{u) = u. Now, if we were to note 
that in the holomorphic domain the R transformation for the Gaussian anti-hermitian 
ensemble is Rmiu) = —u, we readf] 

RiH+hi{u) = R iH {u) + R M (u) = -u + R M {u). (40) 

These two equations yield 

Bi H+ M{u) = B H+M {u) - 2u (41) 

*Note that anti-hermitian Gaussian nullifies in the holomorphic domain the hermitian Gaussian in the 
sense of the group property of additive transformation R, i.e. Rh + Rw = 0. 
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where we have used the relation B(u) = R(u) + 1/u. Substituting u — > Gh+m{z) we 
can rewrite (HIT) as 



B iH+M [G H+M {z)} = z - 2G H+M {z) . (42) 
Let w be a point in the complex plane for which Gih+m{w) = Gh+m(z). Then 

w = z- 2G H+M (z) . (43) 



Equation fl4"3|) provides a conformal transformation mapping the holomorphic domain 
of the ensemble H + M (i.e. the complex plane z minus cuts) onto the holomorphic 
domains of the ensemble iH+M, i.e. the complex plane w minus the "islands" , defining 
in this way the support of the eigenvalues. 

Examples 

Consider the case of "summing" two Hermitian random Gaussian ensembles, i.e. 
consider the Hamiltonian H = Hh + gHn, where g is some arbitrary coupling. The 
sum constitutes of course the Gaussian ensemble, and the spectrum follows from the 
properties of the R function Rn+gni 2 ) — (1 + g 2 )z, or equivalently, Green's function 



G{z) = ^TT?) [z " v' 2 - 4(1 + g2)] (44) 

i.e. the support of the eigenvalues forms the interval (cut) X = [—2^1 + g 2 , +2\/l + g' z \. 
According to ((E|), we can map the interval 1 onto the boundary delimiting the holo- 
morphic domain of the non-hermitian ensemble 7i = Hh + i'yHn, {g — > ij), that 
is 



w = —^\g l z + yjz* - 4(1 + g 2 )} (45) 

i + g v 



with z = t ± iO and t in I. Equation (45) spans an ellipsis with axes 

2/ VI + g' 2 and 



2g 2 / \/l + g 2 . For g 2 = 1 the ellipsis is just the Ginibre's circle. 



A similar construction and an identical mapping (f43D gives the support of the 
eigenvalues in the case of a schematic chiral Dirac operator with chemical potential. Let 
us first consider the case when \i = ie, i.e. the case when the ensemble is hermitian. In 
this case, the resulting Green's function is known to fulfill the so-called Pastur equation 
(random gaussian plus deterministic hermitian Hamiltonian E, here with N/2 levels e 
and N/2 levels -e) 

G ( z ) = 7TT-, + 77T-, ( 46 ) 

v ' z- G(z) - e z - G(z) + e v ' 

encountered in many areas of physics 19 . This is exactly equation (|29|) with the formal 
replacement fi 2 — ► — e 2 . At a particular value of the deterministic parameter e = 1, 
the single cut supporting the spectrum of the hermitian ensemble splits into two-arc 
support, manifesting therefore a structural change in the spectral properties, hence a 
"phase transition" . The spectral properties of the non-hermitian model, with chemical 



potential fj,, follow from the mapping (j43|) , but with Gh+m replaced by an appropriate 



branch of the cubic Pastur equation (46). In particular, at the value fi 2 = 1 the 



spectrum demonstrates the structural change - an island splits into two disconnected 
mirror islands (see Fig.^j). 
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Figure 6. Conformal mappings for the case of Ginibre-Girko ensemble with r = 0.5 (upper), 
non-hermitian chiral ensemble with chemical potential /i 2 = 0.8 (middle) and /z 2 = 1.2 (lower). The 
shaded regions represent the holomorphic domains. 



TWO-POINT FUNCTIONS 

To probe the character of the correlations between the eigenvalues of non-hermitian 
random matrices, either on their holomorphic or nonholomorphic supports, it is rele- 
vant to investigate two-point functions. For example, a measure of the breaking of 
holomorphic symmetry in the eigenvalue distribution is given by the connected two- 
point function or correlator 



N 2 G c (z,z) = ( 



tr 



1 



2 



>c (47) 



z-M 

where the z and z content of the averaging is probed simultaneously. The correlation 



function (47) will be shown to diverge precisely on the nonholomorphic support of the 



eigenvalue distribution, indicating an accumulation in the eigenvalue density. In the 



conventional language of "quarks" and "gluons" , (§7]) is just the correlation function 
between "quarks" and their "conjugates" . A divergence in (|l"7 ) in the z-plane reflects 
large fluctuations between the eigenvalues of the non-hermitian operators on finite z- 
supports, hence their "condensation". 

It was shown in 20 and 21 that for hermitian matrices (with z — > w) the fluctuations 
in connected and smoothened two-point functions satisfy the general lore of macroscopic 
universality This means that all smoothened correlation functions are universal and 
could be classified by the support of the spectral densities, independently of the specifics 
of the random ensemble and genera in the topological expansion (see 22 for a recent 
discussion) . 

In the case on NHRMM the generalized two-point correlator reads 7 

N 2 G c (z, w) = -d z d w logdet(l - g l <g> glT) . (48) 



Here the logarithm is understood as a power series expansion. Equation (|i8 ) is valid for 
Gaussian ensembles and, in the general case, up to factorizable corrections in the sense 
of 10 . The operator Qi (g) £/|T is a tensor product of 2 x 2 matrices (see Fig. |7]). The 
kernel V includes the details of the "gluonic" interactions, depending on the particular 
measure. The tensor structure reflects the nonholomorphic solutions. The choice of 
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fl = £. fl / 2 (8) Sj 152 r/2/355253 

-53 



Si -@- 



Figure 7. Two-point kernel with f,g = q, q. 

"isospin" in the e.g. lower fermion line is done by choosing the appropriate derivative 
d w for the "quark" and c\j for the "conjugate-quark" . 

For the holomorphic solutions of NHRMM, the structure is simpler, since the 
Green's functions are holomorphic in this case. 

Examples 

In the case of the Ginibre Girko ensemble, the two-point correlator in the holo- 
morphic domain reads simply 

C(z, z) = -d z d- z ln[l - G{z)G{z)\ . (49) 

Indeed, in the holomorphic domain instead of the matrix valued Q we use G(z) given 
by (^1|), and the kernel T reduces to the "quark-conjugate quark" coupling equal to 1 
(upper-right corner of the last matrix in (|15|)). Note that the correlator ([49|) diverges 
on the line 

1-|G(^| 2 = (50) 

therefore the ellipse (0), confirming our statement. 

For the nonholomorphic domain, the calculation is a bit more involved, due to the 
explicit matrix structure of Q and T. The explicit form of the matrix- valued resolvent 
is 



1 ( z-rz g z 



1-t 2 \ g z Z-TZ 



(51) 



with g 2 = \z\ 2 (l — r) 2 — r{z + z) 2 — (1 — r 2 ) 2 . One recognizes (^) as the upper left 
corner of (|5T|). The explicit form of the kernel is 

r = diag(r,l,l,r) (52) 

which corresponds to all possible contributions from four graphs on Fig. After some 
algebra, the determinant of 4 by 4 matrix (1 — Qi eg) Q^) turns out to be equal to 
\z — w\ 2 , giving the correlator 



N 2 G qq (z,w) = — 1 (53) 
(w — z) 



In the non-hermitian chiral case the correlator in the outside (holomorphic) region is 
calculated using the same arguments as above. The only minor technical complication 
stems from the chiral (block) nature of the Green's functions. We skip the details 
published elsewhere 7> 15 , showing only the final result. The determinant of 1 — Q\ ® Q\Y 
gives 

(D-fi 2 ) 2 - \z-G\ 4 



D 2 

where the holomorphic G is the appropriate branch of ( p9|) and D = \(z — G)/G\ 2 . 



(54) 
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The zero of the determinant in (pESP occurs for (D — /i 2 ) — \z — G\ 2 , that is 

|^ - C| 2 (l - |C| 2 ) - ^ 2 |C| 2 = (55) 

as quoted in 15 . This is exactly the equation of the boundary separating the holomorphic 
and nonholomorphic solutions, obtained in examples before either as £2 = or as a 
result of conformal mapping. 

In the hermitian case fi = (and z — w), the determinant in ([54]) is simply 
(1 — G 2 (z)G 2 {w)) (chiral) as opposed to (1 — G(z)G{w)) (non-chiral). As a result, for 
w = z and /i = 0, ([54]) is 

N 2 G(z,z) = z 2 J_ A? (56) 

which coincides with (5.5) in 23 . 

Note that the form of fl5"6|) signals two kinds of microscopic universalities. The 
1/N expansion breaks down at z = ±2 (endpoints of spectra) and z = ( "Goldstone" 
pole due to the "chiral" nature of the ensemble). 

The divergence at z = ±2 points at the edge universal behavior of the spectral 
function (Airy oscillations) 24 , the divergence at z = signals the chiral microscopic 
universality 25 . Unfolding the spectra at these points allows to get the explicit universal 
kernels characterizing the fore mentioned universalities. 

In the light of the above remarks, it is tempting to speculate 7 , that the divergences 
of generalized correlators may onset some new types of microscopic universalities present 
in the NHRMM. 

Note also that the relations fl50[), ([54]) demonstrating the functional dependence 
of the two-point correlator on one point holomorphic Green's function allow for some 
extension of the macroscopic universality for NHRMM as well. The eventual geometric 
interpretation of this extension remains an open problem. 

Before closing this example let us present for completeness the result of the chiral 
correlator in the nonholomorphic domain. The calculation is tedious, due to the fact 
that in the nonholomorphic region "quarks" may turn to "conjugate-quarks" and vice- 
versa, with all "quark" species interacting with themselves, and species appear in chiral 
copies. Nevertheless, the final result for the determinant is remarkably simple: 

, , /, rl K~ rTf 1 12, , | 2 (m 2 — -y 2 ) (v 2 - v 2 )) 2 -v 2 y 2 . . 
det(l— y(z) ® y (w)F) = \z—w\ \z+w\ (57) 

(where y = lmz,v = Imw), suggesting perhaps the possibility of further technical 
developments in the case of NHRMM. 



PARTITION FUNCTIONS 

We show now that the information carried by the one- and two-point functions 
is sufficient to specify the "thermodynamical" potential to order 0(1/ N) in the entire 
z-plane modulo isolated singularities, as we now discuss. Similar ideas were used in 
different context in 26 ' 27 ' 28 . 

Let Zn be the partition function in the presence of an external parameter z. In 
the 1/N approximation, the diagrammatic contributions to the partition function Zn 
read 

log Z N = NE + E, + 0(h (58) 
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where Eq is the contribution of the "quark" or "conjugate quark" loop in the planar 
approximation, and E\ is contribution of the "quark-quark" loop, and so on, in the 
same approximation. We will restrict our attention to non-hermitian matrices with 
unitary randomness, in which case the non-planar corrections to Eq are of order 1/N 2 . 
Hence, E is determined by the one-point function and E\ by the two-point functions. 

For z such that (|58|) is real, continuous and nondecreasing function of the extensive 
parameters 29 , log Z^/N may be identified with the "pressure" of the random matrix 
model. As a result, the isolated singularities in the "pressure" are just the "phase" 
boundaries provided that the expansion is uniform. Below we give examples where the 
"phase" boundary is either mean-field-driven or fluctuation-driven. We have to distin- 
guish two cases: holomorphic partition functions ( "unquenched" ) and nonholomorphic 
partition functions, where the complex phases of the determinants are neglected. 

Holomorphic Z 

We consider the partition function 

Z N = (det(> -M)) = Q dtfjchpe-^ {z - M) ^ . (59) 

In contrast to the one- and two-point correlators discussed above, the determinant in 
( |59"D is not singular in the z-plane configuration by configuration. Hence, ( |59"D is a priori 
holomorphic in z (minus isolated singularities). 

The one- and two-point functions on their holomorphic support may be obtained 
from log Zn by differentiation with respect to z. Therefore, from (|5S|) 

E = J Z qIz G(z') + const (60) 

or equivalently 

E = zG- JdG z(G) + const (61) 

after integrating by part. Note that z{G) = B(G) is just the Blue's function 17 of 
G. The constant in E is fixed by the asymptotic behavior of ([59"D, that is Z N ~ z N . 
The planar contribution to Ei in flolf ) follows from the "quark-quark" wheel (two-point 
correlator). The final result for Zn is 

Z N = e NE ° • ({det(l -G®g T -T)Y h + (1)) . (62) 

Note that due to the power —1/2, the fluctuations have "bosonic" character, and are 
dwarfed by the "quark" contribution as (1 : N) 30 . Both Eq and E\ are simple func- 
tions of the resolvent on a specific branch, as expected from generalized macroscopic 
universality. 

We note that the partition function Zn through (|62|) exhibits an essential singular- 
ity in 1/N as expected from thermodynamical arguments. Assuming that the expansion 
for log Zn/N is uniform, then log Z^/N follows from (|6"2]) using the holomorphic resol- 
vent G(z) for large z. The small z region follows by analytical continuation. However, 
since G(z) is multi- valued (already the simple case of the Ginibre-Girko ensemble yields 
two branches for the resolvent in , the analytical continuation is ambiguous. The 
ambiguity may be removed by identifying log \ Z^\/N with some generalized "pressure" 
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and taking G(z) so that log|Zjv|/iV is maximum. As a result, Vn = log Z^/N is 
piece- wise analytic in leading order in 1/N with "cusps" at 

F m ( X) y ) = V j{) ( x? y ) _ V <i) (x, y)=0, (63) 

following the transition from branch % to branch j of G. 

The character of the transition in the 1/N approximation can be highlighted by 
noting that for any finite N, the partition function (|59"D is a complex polynomial in z 
of degree N with random coefficients. In large N, 

1 lr 

V N = — log \Z N \ = - J dvdvg(v,v) log \z - v | 2 . (64) 

To leading order, the distribution of singularities along the "cusps" fl5"3] ) is 

g(z,z) = ±-\d z F\ 2 6(F(z,z)) (65) 

which is normalized to 1 in the z-plane. Redefining the density of singularities by unit 
length along the curve F(z, z) = 0, we may rewrite ( ^5]) as 

Q(z,z)\ F=0 = ±\d z F\ = i-| G » _ G W|. (66) 



For q ^ 0, the integrand in (|64|) is singular at z — v which results into different forms 
for Vn, hence a cusp. For g = 0, that is d z F = 0, Vjv is different iable. For physical 
Vn (real and monotonically increasing), the points g = are multi-critical points. At 
these points all n-point (n > 2) functions diverge. This observation also holds for Ising 
models with complex external parameters 31 . Assuming macroscopic universality 20 
for all rz-points (n > 2), we conclude that d z F = means a branching point for the 
resolvents, hence d z G = oo or B'(G) = 17 . For hermitian matrices, these conditions 
coincide with the end-points of the eigenvalue distributions 17 ' 32 . 

Examples 

To illustrate the above concepts, consider again first the Ginibre-Girko ensemble. 
The resolvent in the holomorphic region satisfies (2~ID, so 

z = rG + i. (67) 



The integration (|6lD in E is straightforward, and after fixing the asymptotic behavior 
we obtain 



1 N = G~ N e^ NG2 ( (1 - G'\z)r)^ + Q I J_ ) ) . (68) 



Here G is the solution of (37). The pre-exponent in (]68|) follows from ( |6^) with the 
matrix Q replaced by G and T = r, as seen in the "quark-quark" component of the 
vertex matrix in (|i~5|). Using ( 67]) we observe that the pre-exponent diverges at two 
points in the z-plane, z 2 = 4r. At these points there is a "phase" change as we now 
show. 

Given (|68D , the generalized "pressure" in leading order is 

V± = -\log (G±G ± ) + T -{Gl + Gl) +o(jjV (69) 
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V± define two intersecting surfaces valued in the z-plane, for two branches G± of the 
solutions to fl67|) . The parametric equation for the intersecting curve is 

F(z, z) = V + -V_ = 0. (70) 

As indicated above, Vjv is piece-wise differentiable. Note that F = on the cut along 
the real axis, —2^/r < z < +2y/r, and from ( |BE| ) the density of singularities per unit 
length is 

1 V4r - z 2 

Q{z,z)\ F=0 = - ^ ' ( 71 > 

Along F, the density of singularities is semi-circle. The density ( [7TD vanishes at the end- 
points z = ±2y/r. This is easily seen to be the same as d z G = oo, or dB{G)/dG = 
with B(G) = tG+1/G. As noted above the term in bracket in Eq. (^) vanishes at these 
points, with a diverging "quark-quark" contribution. The transition is fluctuation- 
driven. These points may again signal the onset of scaling regions with possible uni- 
versal microscopic behavior for non-hermitian random matrix models. This issue will 
be pursued elsewhere. At these points the 1/N expansion we have used breaks down. 

Let us move now to the chiral non-hermitian ensemble. Elementary integration 
for this case leads to 

Z N (z, /x) = e"*» ■ {{D-\D + » 2 ) 2 -(z- G) 4 ]} ^ + O (A)) , (72) 
but now D = (z - G) 2 /G 2 , and 

E (zX) = G 2 + log^- (73) 

with the appropriate branch of holomorphic G solution to (p9|). 

Note that for z = and G 2 = — 1—^j 2 , the pre-exponent in (|72|) diverges. It 
also diverges at z = z* which are the zeros of (|66D for the present case (two zeros 
for small fi and four zeros for large /i), see Fig. ||. Again, at these points, the 1/N 
expansion breaks down marking the onset of scaling regions and the possibility of 
microscopic universality. The z = divergence is just the notorious "Goldstone" mode 
in chiral models, illustrating the noncommutativity of iV — > oo and z — > 0. The rest 
of the arguments follow easily from the preceding example, in agreement with the 
"thermodynamics" discussed in 30 . The analytical results for the nature and location 
of singularities of this example were confirmed by an extensive numerical analysis of 
Yang-Lee zeroes (up to 500 digits accuracy) in 33 . 

Nonholomorphic Z 

The above analysis for the holomorphic thermodynamical potential may also be 
extended to nonholomorphic partition functions of the type 

Z N [z, z] = (det \z -M\ 2 ) = Q d$d$d4dfc^*- M W-+&- M ^*\ . (74) 

Note the important "quenching" of the phase of the determinant in comparison to the 
holomorphic case ( [j)9|) . As a consequence, the two-point correlators diverge rather on 
the one-dimensional boundary separating the phases (when approaching the boundary 
from the holomorphic domain) then at discrete points. 

Similar reasoning as before leads to the explicit expression 

Z N [z,z}=e NEo ■ ({det(l-0 (8) ~f ■ Y) |det(l-£ <g> Q T ■ r)|}"+oQ) (75) 
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where Eq comes from the solutions of 

G(z, z) = G(z, z) = . (76) 

Note that again the contributions from the two-wheel diagrams are of the form 1 / y/ det 
and hence "bosonic" in character. The result could be easily guessed without perform- 
ing the calculations: there are two contributions from the "quark-conjugate-quark" 
wheels (correlators) (square of the 1 / y/ det in first term in the curly bracket) one con- 
tribution from the "quark-quark" wheel and one contribution from the "conjugate- 
quark-conjugate-quark" wheel (represented together as a second (modulus) term in the 
curly bracket). 

Again, the partition function Z N has an essential singularity in l/N, but log Z N /N 
does not. For any finite N, the latter diverges for 

det(l-£®cf -r) = (77) 

which defines the line of singularities, and for 

det(l - Q ® Q T ■ T) = (78) 

defining the set of discrete points, encountered in the case of the holomorphic partition 
function. 

Examples 

For the Ginibre-Girko example, the line of singularities ( |77"D reads 

1-\G\ 2 = 0. (79) 

The line of singularities (|7|) reproduces in this case the ellipse (p3|). The ellipse includes 
the points of the "phase" change (see Q5"8|)), 

l-rG 2 (^) = (80) 

i.e. the focal points z 2 = 4r, corresponding to (|78|), connected by the interval (|70|) , i.e. 
F = 0. 

In the case of chiral non-hermitian random model the condition (|77|) reads 15 

D~ 2 [(D - jj 2 ) 2 -\z- G| 4 ] = (81) 

with D = \(z — G)/G\ 2 , therefore exactly the condition (|55|). This line represents the 
boundary between the holomorphic and nonholomorphic solutions. The set of discrete 
multi-critical points, corresponding to flTBD is given by the condition 



D- 2 [(D + fi 2 ) 2 - (z - G) 4 ] = (82) 

but with D = [(z — G)/G] 2 , in agreement with fl7"2]). Note the crucial appearance of the 
modulus and the flip in the sign of fi when comparing last two formulae. The explicit 
solution of ( |52"[ ) consists on set of two or four points, (depending on the value of the /j,), 
being the analogs of Airy type end-points singularities and a single multi-critical point 
z = 0, reflecting the chiral nature of the ensemble. 

Figure [8] shows the critical lines and critical points corresponding to the conditions 
(|77| , |78D for Ginibre-Girko and non-hermitian chiral ensembles. End-points singularities 
are denoted by "NATO stars" , chiral singularity - by "Zakopane sun" . 
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Figure 8. Critical lines (f77j) and critical points (|7S| ) for Ginibre-Girko ensemble with r = 0.1 (left) 
and non-hermitian chiral ensemble with small (/i 2 = 0.1, middle) and large (/i 2 = 0.6, right) value of 
chemical potential. The solid lines and the symbols represent the manifold on which the 1 /N 
expansion breaks down and signal the location of possible microscopic universal behavior for 
NHRMM. 



The fact that the critical line in Fig. ^ b,c surrounds the multi-critical points of 
the unquenched partition function, explains the failure of quenched lattice calculations 
with finite baryonic potential. The nature of chiral restoration is masked by unphysical 
fluctuations caused by neglecting the phase of the determinant. The critical line (|8lD 
exactly reproduces the shoreline of islands of unphysical "mixed-condensate" , obtained 
using the replica methods 9 . 



CONCLUSIONS 

Most of the details of results presented here are included in already published 
papers 6 ' 7 ' 15 . In this mini-review on some novel aspects of NHRMM, we tried to 
enhance the role of a priori not obvious connections between the hermitian and non- 
hermitian ensembles of random matrices. In particular, we presented several ways 
for providing the supports and the density of eigenvalues for non-hermitian ensembles 
and the way for calculating smoothened (wide) correlations, via either analogies (di- 
agrammatic expansion, generalized Blue's functions), or formal relations between the 
hermitian and non-hermitian ensembles (conformal mapping, quenching/unquenching 
of partition function). We used the same set of known examples to demonstrate clearly 
the cross-references between the methods, as well to exhibit the shortcomings of the 
standard treatment. 

We also speculated on some new features, like generalization of macroscopic uni- 
versality and the possibility of several types of new microscopic universalities related 
to the critical behavior of various correlators in the NHRMM. 

This last issue is of great interest in light of recent exciting results in "weakly" 
non-hermitian random matrix models 34 and related applications 35 . 
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